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We consider orthogonal, unitary, and symplectic ensembles of random matrices with (l/a)(lnx) 2 
potentials, which obey spectral statistics different from the Wigner-Dyson and are argued to have 
multifractal eigenstates. If the coefficient a is small, spectral correlations in the bulk are universally 
governed by a translationally invariant, one-parameter generalization of the sine kernel. We provide 
analytic expressions for the level spacing distribution functions of this kernel, which are hybrids of 
the Wigner-Dyson and Poisson distributions. By tuning the single parameter, our results can be 
excellently fitted to the numerical data for three symmetry classes of the three-dimensional Anderson 
Hamiltonians at the metal-insulator transition, previously measured by several groups using exact 
diagonalization. 
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I. INTRODUCTION 

Quantum mechanics described by stochastic ensembles 
of Hamiltonians |l| , and by Hamiltonians with classically 
chaotic trajectories 0, have been a subject of intense 
study for years. In contrast to Hamiltonians of classi- 
cally integrable systems whose energy levels are mutu- 
ally uncorrelated, chaotic Hamiltonians generally exhibit 
strong correlation among levels. To simulate this level 
repulsion in chaotic or disordered Hamiltonians, Wigner 
introduced the random matrix ensembles (RME) In 
the RME defined as an integral over N x N matrices, 
only the antiunitary symmetry of the Hamiltonian is re- 
spected, and its spatial structure is completely discarded. 
Despite this extreme idealization, analytic predictions 
from RME's beautifully explain the spectral statistics of 
chaotic or disordered Hamiltonians Q. This success is 
accounted for by the fact that the RME is not merely an 
idealization of disordered Hamiltonians but it is indeed 
equivalent to the latter under a situation where the mean 
energy level spacing A is much smaller than the Thouless 
energy E c (inverse classical diffusion time) that the di- 
mensionality does become unimportant This on the 
other hand implies that in a region where the mean level 
spacing is equal to or larger than the Thouless energy, 
and the associated states tends to localize due to diffu- 
sion, standard RME's cannot provide good quantitative 
descriptions. 

As a concrete example of systems of disordered conduc- 
tors, le_t us take the Anderson tight- binding Hamiltonian 
(AH) ■ 



e r aia T 



<r,r') 



(1) 



describing free electrons in a random potential. Here 
a\ and a r are creation and annihilation operators of an 
electron at a site r on a three-dimensional (3D) toroidal 
lattice of size L 3 , e r 's are site energies that are mutually 



independent stochastic variables derived from a uniform 
distribution on an interval [— W/2, W/2], and the sum 
(r,r') is over all pairs of nearest-neighboring sites. The 
antiunitary symmetry of the AH [Eq. (Q)] is orthogonal 
(/3 = 1), as it respects the time-reversal symmetry and 
has no spin dependence. For small values of disorder W, 
all eigenstates are extended. An overlap of such extended 
states \l/i(r) and \&j(r) is expected to induce repulsion be- 
tween associated eigenvalues, of the form \xi — Xj\. Then 
the statistical fluctuation of these eigenvalues is described 
well by the Gaussian orthogonal ensemble of random ma- 
trices whose joint probability distribution consists of the 
product of a Vandermonde determinant rii<j \ x i ~ x j\ 
and Gaussian factors Yii e ~~ Xi ■ The quadratic potential 
in the latter is merely for the sake of technical simplicity, 
and the deformation of the potential by generic polyno- 
mials universally leads to the same Wigner-Dyson statis- 
tics ||. As we increase the disorder W, states associated 
with eigenvalues close to the edges of the spectral band 
are believed to start localizing. When the disorder is as 
large as to induce the metal-insulator transition (MIT) 
E c / A ~ 1, eigenstates are observed to be multifractal 
Q, characterized by an anomalous scaling behavior of 
the moments of inverse participation ratio jio] |l2| 



(2) 



This property implies a slowly decreasing overlap be- 
tween states |T|] (for \xi — Xj\ ^> A), 



(3) 



For large enough disorder, the off-diagonal (hopping) 
term of the AH becomes negligible, leading to mutually 
uncorrelated eigenvalues. Each cigenstate is almost local- 
ized to a single site. Thus the spectrum of the AH shows a 
gradual crossover from the Wigner-Dyson to the Poisson 
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as one increases the disorder W, while keeping the size 
L fixed finite. A characteristic observable in these stud- 
ies of spectral statistics is the probability E(s) of having 
no eigenvalue in an interval of width s, or equivalently 
the probability distribution of spacings of two consecutive 
eigenvalues P(s) = E"{s). These observables capture the 
behavior of local correlations of a large number of energy 
eigenvalues, as the former consists of an infinite sum of 
integrals of regulated spectral correlators (the subscript 
reg denotes its regular part, i.e., with ^-functional peaks 
at coincident x,'s subtracted), 

~ f s / 2 
E ( S ) = Z2^^~ / dx 1 ---dx n (p(x 1 )---p(x n )) tcg ,(4 : ) 

and are more conveniently measured by the exact di- 
agonalization of random Hamiltonians than the spec- 
tral correlators. The Poisson distribution is character- 
ized by -Pp(s) = exp(— s), and the Wigner-Dyson dis- 
tribution is well approximated by the Wigner surmise 
P w (s) = (vrs/2)exp(-7rs 2 /4)- 

Recent technical developments [[l4]-p6[ on the exact 
diagonalization of the AH on a large size of lattices 
prompted analytical studies on its spectrum J27|-|3l|]. It 
was noticed in Ref. jL4j that at the MIT point with dis- 
order W ~ 16.5, the level spacing distribution function 
(LSDF) P(s) is independent of the size L. From this 
finding these authors have argued that in the thermo- 
dynamic limit there exist only three universality classes: 
Wigner-Dyson, Poisson, and the third, critical statistics. 
The presence of critical LSDF's was also observed for 
the unitary (j3 = 2) and symplectic (8 — 4) cousins, i.e., 
AH's under a magnetic field , 

H = 'S^e r ala T + V rr >ala T >, (5) 

r <r,r') 
V ,~ =e T2Kiar y y = y _ _ 1 

and with spin-orbit coupling |I(|,|25|,^lJ , 

H = e t al a a ra + ^ V Tafr ' a >al a a T * a >, (6) 

r,a=± (r,r'),CT,cr' 

V Ta , r± u, = (e^ 1 )^ (i = x,y,z). 

The critical LSDF's are found to be independent of the 
strength of the magnetic field a or the spin-orbit cou- 
pling 9. For all values of /3, the Wigner-Dyson-like be- 
haviors P< CT ) (s) oc s@ for small s have been confirmed. 
There were disputes over the large s asymptotic behav- 
ior of the LSDF's, but accurate measurements on large 
lattices (2^] strongly support the Poisson-like behavior 
P( cr -j(s) ~ exp(— const x s) for large s, excluding a non- 
trivial exponent P( cr )(s) ~ exp(— const x s 1+7 ) predicted 
in Ref. p9] . In the light of the success of the random ma- 
trix (RM) description of extended states, a natural resort 
to describe this critical statistics is to consider a deformed 
RME that violates the above mentioned universality of 



the Wigner-Dyson statistics for the Gaussian ensembles. 
The validity of the RM description for such critical sys- 
tems is of course far from clear, because the existence of 
the MIT crucially relics upon the dimensionality, whereas 
the RME has no spatial structure. Nevertheless, under 
an assumption that the spectra of the AH be described 
by RME's, an attempt |32| was made to reconstruct a 
random matrix potential out of the macroscopic spectra 
of the AH. There, it was observed that a potential of the 
form 

V{x) = UHl + b\x\)] 2 (7) 
la 

explains well the numerical data. The above potential in- 
deed violates the universality of the Wigner-Dyson statis- 
tics that is guaranteed only for polynomial potentials 
p3[ . This observation leads to a speculation that the 
critical level spacing distribution might be derived from 
a RME with a potential (0). Later the LSDF of orthog- 
onal RME's of type (0) has been measured by using the 
Monte Carlo simulation |?4| in order to compare it with 
that from exact diagonalization of the AH p2|. Excel- 
lent agreement between the two was found by tuning the 
parameter a to 2.5. Motivated by this success, we shall 
derive analytic forms of the LSDF's of RME's with or- 
thogonal, unitary and symplectic symmetries. (The uni- 
tary case has already been reported |3f|.) In doing so, 
we shall retain all perturbative (polynomial in a) parts 
of the spectral kernel, and discard unphysical nontrans- 
lationally invariant parts of order 0(e _7r / a ). 

This paper is organized as follows. In Sec. II we review 
nonstandard features of RME's with log-squared poten- 
tials. In Sec. Ill we follow the method of Tracy and 
Widom [^7j to derive the LSDF's from an approximated 
translationally invariant kernel. In Sec. IV we shall com- 
pare our results with the numerical data of the AH's with 
orthogonal, unitary, and symplectic symmetries. In Ap- 
pendix A we collect standard results on Fredholm deter- 
minants in random matrix theories that are relevant for 
our purpose. 

II. RME WITH LOG-SQUARED POTENTIAL 

In this section we review properties of RME's with 
log-squared potentials. For small enough a, we derive a 
translationally invariant kernel, whose level spacing dis- 
tribution will be our subject in this paper. 

We consider TV x N random real symmetric (f3 = 1), 
complex hermitian (/3 = 2), and quaternion selfdual 
(/3 = 4) matrix ensembles, whose joint probability densi- 
ties of eigenvalues are given by 

N 

^(A 1 ,...,A A r)cxne- 1/(A * ) I] \ X i- X jf- (8) 

i=l l<i<j<N 

with a potential growing as 



2 



y(A)~-(lnA) 2 (A»l). 
la 



(9) 



The potential is assumed to be regularized at the origin, 
as in Eq. ((?]) . For a particular form of the potential 



V(X) 



ln[l + 2<f cosh(2 arcsinhA) + q 2n ] , (10) 



(0 < q < 1), which behaves as Eq. (^) with a = ln(l/g)/2, 
corresponding orthogonal polynomials are known as the 
q-Hermite polynomials [j36| . Using their asymptotic form, 
Muttalib et al. p7j] have obtained the exact kernel for the 
unitary ensemble with the potential ( |lo|) in the large N 
limit, 



K^(x,y) = const x V cos ^ cos ^ x 

cosha(x + y) 

tf4{x + y,e-* 2 / a ) d 1 (x-y,e- 7 ' 2 / a ) 



^di(2x, e-* 2 / a )-di(2y, e-* 2 / a ) sinha(x-y) 

Here "&„ are the elliptic theta functions, and the variable 
x = A/ (2a) is so rescaled that the average level spacing is 
unity. In order to eliminate nonuniversal effects involved 
by the regularization of the potential in the vicinity of 
the origin from the bulk correlation, we need to take 

x, y 3> 1, with | a; — y\ = bounded. (12) 

Then we obtain an asymptotic form of the kernel 

K^ sympt \x,y) = const x 

M x + y,e~* 2/a ) ^i(x-y,e" 7r2 / a ) 



a/i9 4 (2x, e~ n2 / a )i} 4 (2y, e~ 7 ' 2 / a ) sinh a(x - y) 



(13) 



The above kernel is still not translationally invariant due 
to the reason explained below. Now we make a further 
simplification of the kernel by using an approximation. 
For e-* 2 / a < 1, we can discard subleading orders from 
the q expansion of the theta functions in terms of trigono- 
metric functions. Then we obtain a translationally invari- 
ant kernel, up to 0(e~ 7T / a ), 



K(x,y) 



a sin7r(a; — y) 
7r sinha(izi — y) 



(14) 



The universality of this deformed kernel within RME's 
is observed for the g-Laguerre unitary ensemble [j38| , 
the finite-temperature Fermi gas model |3j|, and subse- 
quently for unitary ensembles whose potentials have the 
asymptotics ([)]) |Q . This universality can be considered 
as an extension of the universality of the sine kernel H , 



K(x,y) = 



sin7r(a; — y) 
ir(x-y) 



(15) 



for polynomially increasing potentials, proven via the 
asymptotic WKB form of the wave functions 



tp N (X) ~ cos ^tt J p(X)d\ + , 
K( A ,A0~- n[7r(/Ap ~ /A ' P)] 



A - A' 



(16) 
(17) 



Here p(X) stands for the exact unnormalized spectral den- 
sity, if (A, A). In the spectral bulk of the RME with 
polynomially increasing potentials, the spectral density 
divided by N is bounded, and is locally approximated by 
a constant when measuring A in unit of the mean level 
spacing. This slowly- varying function is called the mean- 
field spectral density p(A), given by GO] 



N 



1 f R dp R 2 -X 2 V'{p) 



ttVR 2 - A 2 ^ 2 J-R X - /i I/ R 2 - M 2 2 



(18) 



Here ±i? are the end points of the spectrum, determined 
by J R R p(X)dX = N. After replacing the exact p(X) by 
the mean-field p(A), the unfolding map 



A i— > x = / p(X)dX 



(19) 



becomes merely a linear transformation, leading univer- 
sally to the sine kernel. On the other hand, in our case 
of the potential ([}]) , the mean-field spectral density jl^ ) 
behaves as 



i 



2a|A|' 



implying an unusual unfolding map 



1 



A i— > x — — sgn(A) In |A| 
2a 



(20) 



(21) 



while the formula (|17|) stays valid |4l). Then the kernel 
( |l7| ) universally reduces to Eq. ( |l4| ) after this unfolding. 

It is clear from the form of the kernel (14) that a set 
of eigenvalues with \xi — Xj\ 3> 1/a obeys the Poisson 
statistics, i.e., is uncorrelated. On the other hand, a set 
of eigenvalues with \xi — Xj\ <C 1/a obeys the Wigner- 
Dyson statistics, because Eq. ( [fi]) is then approximated 
by the sine kernel, up to 0(a 2 ). To be precise, the kernel 
(Q) signifies the mult ifr act ality of the eigenstates [EgJ . 
To see this, we first note that the property (|3|) of the frac- 
tal states leads to a compressible gas of eigenvalues, i.e., 
a linear asymptotics of the number variance X 2 within an 
energy window of width S W3[ [^(a^) = 1 — (/°( a; )/ o (0))reg], 



r s 

E 2 (S') = S - 2 / dx(S - x)Y 2 (x) 
Jo 



(22) 



Ul-^)S X S (,s 1). (2.-!) 
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The RME's with the scalar kernel ( |i4| ) indeed enjoy 
this asymptotic behavior with the level compressibility 
X given by 



X 



0(e 



-2ir 2 /a\ 



(24) 



for unitary |38|,[42| and orthogonal ensembles, and for 
the symplectic ensemble the exponential correction is re- 
placed by 0(e _7r '"). Moreover, the mii/tifractality of 
the RME's with the deformed kernel ( |l4|) has been con- 
cluded |Q through the equivalence between the finite- 
temperature Fermi gas model (having the universal de- 
formed kernel) and a Gaussian banded RME that is 
proven to have multifractal eigenstates |44|. This brings 
forth the possibility of describing the spectral statistics 
of the AH's at the MIT by RME's with the deformed 
kernel. 

We should emphasize an important fact that the 
RME's with log-squared potentials cannot describe dis- 
ordered systems with large disorder. While the AH in the 
W — > oo limit leads to the Poisson statistics, the RME 
with Eq. (||) does not obey the Poisson statistics in the 
limit a — > oo pi| , fl5| ] . It is because the joi nt p robability 
distribution of RME's after the unfolding (21) leads to 



Vfj (xi , ■ ■ ■ , xn) = const x 



N 



l<i<j<N i=l 



In the limit a — * oo, each factor of | e 2ct I ' zpe 20 ^^] 
is dominated by an exponential with a larger modulus. 
Thus the Vandermonde determinant is approximated by 



N 



2ali{i-l)\xi 



(for |aci| < •■• < \x N \). 



(26) 



Consequently Vp{x\, . . . , xn) tends to a product of very 
narrow [of variance a 2 = l/(4a/3) -C 1] Gaussian distri- 
butions obeyed by Xi whose center is at [j3(i — 1) + l]/2. 
This "crystallization" of eigenvalues invalidates naively 
expected mutual independence of distributions of eigen- 
values, and drives the spectrum toward an exotic statis- 
tics different from Poissonian pl| . This phenomenon can 
be rephrased in the context of using the WKB formula 
(||) to derive the kernel Although Eq. @ remains 
valid even in the case a — ► oo, use of the mean- field spec- 
tral density p~(X) = l/(2aA) in place of the exact p(X) 
is not justifiable, because the crystallization of eigenval- 
ues leads to a rapidly oscillating p(X). In the case of the 
g-Hermite ensemble ( |l0| ) , the potential itself has a oscil- 
lation of the same type, leading again to crystallization 
pEfl . Therefore, the RME with log-squared potentials, 
despite the fact that it is constructed from the macro- 
scopic spectra of the AH, should be considered as a good 
model of the latter only for small values of a where the 



level repulsion property of the Gaussian ensembles is de- 
formed slightly but not to the extent that the crystal- 
lization of eigenvalues becomes prominent. In the case 
of the q-Hermite ensembles, we can estimate this scale 
to be characterized by the value of a where nontransla- 
tional invariance of the exact kernel becomes manifest, 
i.e., e -7r "/ a ~ 1. We assume this estimate to be valid 
generically for RME's with log-squared potentials (^). In 
view of this, for our purpose of reproducing the LSDF of 
the AH, we shall concentrate on the approximated uni- 
versal kernel (|l4|). This will be justified a posteriori, after 
confirming that the best-fit value of the parameter a for 
the MIT point of the AH's are such that e _7r2/a < 1. 



III. LEVEL SPACING DISTRIBUTIONS OF THE 
DEFORMED KERNEL 

In this section we analytically compute the LSDF 
Pp(s) from the deformed kernel (|1J) for all values of the 
Dyson index /3. Our result completes earlier attempts to 
compute P2(s) numerically (37) or asymptotically [^6| , 
and is consistent with those. 

We notice that the kernel (|l4|) is equivalent to that of 
Dyson's circular unitary ensemble at finite N M, 



K(x,y) = 



sm(N/2)(x-y) 



Nsm(l/2)(x-y) , 
by the following analytic continuation 

2a 



a 



(27) 



(28) 



Tracy and Widom |4j have proven that the diagonal 
resolvent kernel of Eq. (E7j) is determined by a second- 
order differential equation that is reduced to a Painleve 
VI equation |j48|l. We shall follow their method below. 
The kernel (14J) is written as 




<p(x)tp(y) - ip(x)<f>(y) 



(29) 



tp(x) = 



These component functions satisfy 



<fi = a(j> + nip, %l> — —TT(p + aip. 



(30) 



We use the bra-ket notation <f>{x) — (x\(j>) and so forth 
p0| . Due to our choice of the component functions to 
be real valued (unlike jl7|], Sec. V D), we have {x\0\<f>) = 
((f>\0\x) and a similar situation for ib with any self-adjoint 
operator O and real x. Then Eq. ( |29| ) is equivalent to 



[e"*,if] = |0)<V>|-|V>M, 



(31) 



where X and K are the multiplication operator of the 
independent variable and the integral operator with the 
kernel K(x,y)9(y — ti)9(t 2 — y), respectively. Below we 
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will not explicitly write the dependence on the end points 
of the underlying interval [ti, t 2 }. It follows from Eq. (|3l| ) 
that 



,2aX 



K 



'1-K 



1 



1 -K 



- mm) 



1 



1-iT 



(32) 



that is, 



J2ax _ 2ay 



)R(x, y) = Q(x)P{y) - P(x)Q(y), (33) 
Q(x) = {x\(\-K)- 1 \4>), 
P{x) = {x\{l-K)- 1 \i>). 



At a coincident point x = y we have 

2ae 2ax R(x, x) = Q'(x)P(x) - P'{x)Q(x). 
Now, by using the identity 

^ = {-iyK\U){ti\ (i = l,2), 



we obtain 



dQ(x) 
dP{x) 



i-l^Rix, ti)Q(U), 
[-l)*R{x,U)P(U)- 



(35) 

(36a) 
(36b) 



On the other hand, by using the identity (D is the deriva- 
tion operator) 



[D,K\=K(\t 1 )(t 1 \-\ta){t a \), 



(37) 



which follows from the translational invariance of the ker- 
nel (d x + d y )K(x — y) = 0, we also have 



dQ(x) 
dx 



dP(x) 
dx 



--(xlDil-K)- 1 ^) 

-{x\(l-K)-i\ct>') 
+ {x\{l-K)- 1 [D,K]{I-K)- 1 \^) 
-aQ(x) + ttP{x) + R(x, *i)Q(*i) - R[x, t 2 )Q(t 2 ), 

(38a) 

■—irQ(x) + aP(x) + R(x,h)P(h) - R(x,t 2 )P(t 2 ) 

(38b) 



Now we set ti = —t, t<i = t, x, y = —t or t, and introduce 
notations q = Q(—t), q = Q(t),p = P(—t),p = P(t), and 
R = R(-t,t) = R(t,~t),R = R(t,t) = R(-t,-t). The 
last two equalities follow from the evenness of the kernel. 
Then Eqs. (E33|) and (|3^) read, after using Eqs. (|38|), 



pq — qp = 2i?sinh 2at, 
2 ~ 

p 2 + q 2 = - (R 2 sinh 2at + R a e~ 2at ) , 

7T 

2 - -(R 2 sinh2ai + Rae 2at ). 

7T 



2 1 

p +q 



(39a) 
(39b) 

(39c) 



The total t derivatives of Eqs. (p9T) lead to (• = d/dt) 



pp + qq = —(R sinh 2at) ' 
R = 2R 2 , R = ARR. 



(40) 
(41) 



The left-hand sides of Eqs. (|3jj) and (^) satisfy an ad- 
ditional constraint, 

{pp + qq) 2 + (pq - qp) 2 = if + f)(p 2 + q 2 )- (42) 



By eliminating p, p, q, q, R, and R from Eqs. (|39|)-(|f2|) 
we obtain for R(s) (s = 2t,' = d/ds) 



., . sinh as ... . 
a cosh asR'{s) H — R" (s) 



-[tt sinh as R'(s)} 2 = 



(34) R'(s) ([aR(s)] 2 +asmh2asR(s)R' (s) + [sinhas R' (s)} 2 ) . 

(43) 

It is equivalent to Eq. (5.70) of Ref. |47[] after the an- 
alytic continuation (p8|), accompanied by a redefinition 
R(s) — > (i/2a)R(s). This is slightly nontrivial because 



Ref. 1 47 1 has used p — p* and q = q* , which follow 
from the analytic properties of its component functions, 
(f)(—x) = (p(x)* and x) — ip(x)* . In the limit a — ► 0. 
it clearly reduces to the a form of a Painleve V equation 

R'(s) + ^ R "( s )] 2 + ksi?'(s)] 2 = R'(s)[R(s) + sR'(s)} 2 1 

(44) 

derived for the sine kernel ( |l5| ) of the Gaussian ensem- 
bles [Q. Note that for (3 = 1, a replacement a — ► a/2 is 
necessary because of our convention (Al). 

Finally, the LSDF's Ppis) are expr essed i n ter ms of 
the diagonal resolvent R(s) via Eqs. (A12), (A14), and 
the first of Eq. @: 



p i(s) 

P2(S) 
P 4 ( S ) 



e -i f°ds[R( s ) + y/R'(s)] 
- f dsR(s) 

e Jo w 



e a Jo % ' cosh - 
2 



2.s- 



ds^/R'{s) 



(45a) 
(45b) 
. (45c) 



The boundary condition follows from the expansion 
of E 2 (s) in terms of the correlation functions: 



ns/2 

#2(s) = 1- / dxK{x,x) 

J-s/2 
,/2 



dx\dx 2 det K(xi,Xj) 



R(s) 



-s/2 <.3'=1.2 

1 - a + 0(s 4 ), 
-[ln£ 2 (s)]' = 1 + S + -- 



(46) 
(47) 



Our main result consists of Eqs. (]43j), (|45|), and (J47|). 

For s <C 1/a, we can Taylor-expand hyperbolic func- 
tions, to obtain a perturbative solution to Eq. (pl3|), 



5 



R(s 



) = 1+8+8 , + ( 1 _^) 8S + ( 1 _^ ? !)) sl + ( 1 _fel± 



a 2 )(75 - 47r 2 -6a 2 )\ 5 
s - 



450 



Px(s) 

P2(s) 
P 4 (s) 



4tt 2 +cl 2 (4^ 2 +a 2 ) (I2 7r 2 + 7a 2 ) 3 (tt 2 + a 2 ) (4 tt 2 + a 2 ) 



24 



2880 



1080 



(4 tt 2 + a 2 ) (48 tt 4 + 72 tt 2 a 2 + 31 a 4 ) s 5 (^ 2 + a 2 ) (4 tt 2 + a 2 ) (12 tt 2 + 13 a 2 ) 



322560 



226800 



7T 2 + a 2 2 (^ 2 +a 2 )(2^ 2 +3a 2 ) 4 (^ 2 + a 2 )(^ 2 + 2a 2 )(3^ 2 + 5a 2 ) 6 
3 S 45 S + 945 S 

16 (tt 2 + a 2 ) (tt 2 + 4 a 2 ) s 4 128 (tt 2 + a 2 ) (tt 2 + 4 a 2 ) (3 vr 2 + 13 a 2 ) s 6 



135 



14175 



(48) 

(49a) 
(49b) 
(49c) 



The above perturbative expansions are correct for the 
kernel ( p"3| ) of the q-Hermite ensemble in any polynomial 
orders of a, while they lack nonperturbative terms of or- 

2 i 

der e _7r /a , which depend on the reference point. On the 
other hand, for s» 1/a, it can be proven from Eq. ( pl3| ) 
that R(s) approachs a constant. Then Eqs. (Eq) imply 



lnPi(s)~-i J R f (oc) S , 
lnP 2 ,4(s) ~ -P tt (oo)s, 



(50a) 
(50b) 



for s — > oo. In Fig. 1 we exhibit the decay rate n{a) = 
i? a (s = oo) for < a < 4 computed numerically from 
Eq. (^). At present we could not find an analytic form 
of /c(a). For small a(< 0.5), it is well approximated by 
l/n(a) 0.202a, which agrees extremely well with the 
value 2/-7T 2 = 0.2028 expected from Eq. (24) and the an- 
alytic formula [E7| that holds in generality, 



1 

X ~ 2k' 



(51) 



Eqs. (J49j) and (|50|) tells that our LSDF's are indeed hy- 
brids of the Wigner-Dyson-like [Pp(s) ~ s' 3 for s small] 
and the Poisson-like distributions [Pp(s) ~ e _KS for s 
large]. In Figs. 2-4 we exhibit plots of the LSDF's Pp(s) 

2 I 

for p = 1,2,4 and for various a such that e _7r < a -C 1, 
obtained by numerically solving Eq. ([43|). 



IV. ANDERSON HAMILTONIANS AT MIT 

In this section we make comparison between the 
LSDF's and the level number variance in the exact di- 
agonalization of the AH's and our analytic results from 
multifractal RME's. 

As numerical data to compare with, we adopt Ref. f23j | 
for the AH [Eq. (|)] (J3 = 1), Ref. [£§ for the AH un- 
der a magnetic field @ with a = 1/5 (/? = 2), and 
Ref. [g5j for the AH with spin-orbit coupling (||) with 
6 = 7r/6 (fi = 4), at their MIT points. We choose the 
best fit values of a from the exponential decay rates k of 
the numerical data using Fig. 1. Based on the numerical 
results k = 1.9, 1.8, 1.7 for = 1, 2, 4, respectively, we es- 
timate the parameter a in the potentials of RME's to be 



a = 2.95, 3.55, 3.90. Preference could alternately be put 
on best matchings in the smaller values of s (< 2), which 
would lead to a = 3.2 for [3 = 1 [although the difference 
in P\(s) between a = 2.95 and a = 3.2 is tiny]. In Figs. 
5-7 we exhibit linear and logarithmic plots of LSDF's of 
the RME's and the AH's. The numerical data fit excel- 
lently with our analytic result from the kernel (^) , for a 
large energy range < s < 6 where the LSDF's vary by 
four to five orders of magnitude. The use of this approx- 

2 / 

imated kernel is justifiable because e _7r < 1 holds for 
these values of a. Small systematic deviations can be at- 
tributed to the errors involved in determining the values 
of a from the decay rates of numerical LSDF's, and possi- 
bly to an essential difference of order 0(e~™ / a ) between 
the RME's and the AH's. 

Furthermore, we exhibit in Fig. 8 the number variance 
S 2 (S') of the orthogonal AH obtained in Ref. |16|, to- 
gether with the RME result Q with (/? = 1) 



Y 2 (x)=K(x) 2 +K'(x) K(y)dy 



(52) 



at a = 3.2. We can confirm that not only the asymptotic 
slopes x of S 2 (S f ) (first pointed out by Canali |34|| , who 
computed x by the Monte Carlo simulation of RME's), 
but their full functional forms are in a good agreement 
for L < 10. To recapitulate, we have the following three 
distinct functional observables (consisting of the corre- 
lation functions in the second column) that agree well 
between the critical AH's and the deformed RME's: 



Quantity 



Correlation function 



Potential V(A) 
Number variance S 2 (S') 
Level spacing P(s) 



1-level 

unfolded 2-level 
unfolded n- level (n > 2) 



In addition, both the critical AH's M and the deformed 
RME's are shown to have multifractal eigenstates, 
although the sequences of the multifractal dimensions 
are yet to be compared. Agreements in the unfolded 
quantities should not be considered a tautological conse- 
quence of the first line; one should recall that an identical 
semicircle spectrum could as well be obtained either from 
invariant RME's with Gaussian potentials (obeying the 
Wigner-Dyson statistics), or from diagonal random ma- 
trices whose entries are independently derived from the 
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semicircle distribution (obeying the Poisson statistics), or 
from any intermediate ensembles. From these grounds, 
we conclude that the interaction between unfolded en- 
ergy levels of the 3D AH on equilateral (L x = L y = L z ) 
toroidal lattices at the MIT is very well described by the 



form |e 



lax Jlax 1/3 



which is common to the RME's with 



log-squared potentials (||), in contrast to the standard 
form \x — x'\P of the Gaussian RME's and the AH's in 
the metallic regime. We surmise that the dimensionality 
and the fractal dimensionality enter the critical spectral 
statistics primarily through a single parameter a, as long 
as the multifractality of the critical wave functions is not 
too strong. Further work is needed to explain this form 
of the level repulsion from the multifractality (^) of the 
wave functions, and its origin from microscopic models. 

Finally, remarks related to novel numerical results on 
critical AH's are in order. Recently it was observed that 
the critical LSDF of the 3D AH is sensitive to the ge- 
ometr y o f the lattice, i.e., the topology (boundary condi- 
tion) |5l],|5^] and the aspect ratio due to the coher- 
ence of the critical wave functions maintained over the 
whole lattice. Since RME's treat randomness on all cites 
and bonds on equal footing, it is likely that our RME's 
describe best the critical AH's on maximally symmetric 
lattices, i.e., equilateral tori, but not those on less sym- 
metric lattices, such as unequilateral toroidal lattices or 
lattices with boundaries. Besides, since the validity of our 
expressions for the LSDF's is limited to the case of weak 
multifractality (relatively small a), it will not properly 
describe the critical orthogonal AH in four dimensions 
Q, where the level compressibility was observed to be 
larger than the value in three dimensions (x ~ 0-27 |22|) 
and close to its upper bound \ — 0-5- 
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APPENDIX A: FREDHOLM DETERMINANT IN 
RANDOM MATRIX THEORY 

In the Appendix we collect known results on random 
matrix theories that are relevant to our purpose of evalu- 
ating the LSDF of three symmetry classes of RME's. We 
follow Mehta's classical book ||] and the works by Tracy 
and Widom [|^J55|5^]. Readers are referred to them for 



detailed proofs. Subsequently we shall concentrate on 
the case where the spectral correlation is translationally 
invariant after unfolding. 

The joint probability densities of eigenvalues of N x N 
random real symmetric (j3 = 1), complex Hermitian 
((3 = 2), and quaternion self-dual (f3 = 4) matrices are 
given by 

N 

Vp(xi, ■■■ ,x N ) = const x JJw^a;,-) Yl \xi-Xjf, 

i=l l<i<j<AT 

(Al) 



Wl ( x ) = e ~ v( - x)/2 , w 2A {x) = e~ v( - x \ 



We use the above convention between the weight func- 
tions and the potentia ls s o as to simplify the relations 
between kernels [Eq. ( |A5| ) below]. We introduce the 
"wave functions" {V , »(^)}»=o,i,... by orthonormalizing the 
sequence {x l e~ v ^ x ^ 2 }, and the projection operator K 
to the subspace spanned by the first N wave functions. 
As an integration operator acting on the Hilbert space 
spanned by the wave functions, K is associated with the 
kernel (we shall use the same letter for an operator and 
the kernel associated with it), 



JV-l 



K (x,y) = ^2 ipi{x)il>i(y). 



(A2) 



i=0 



Then the joint probability densities are expressed in 
terms of determinants of the kernels 101 : 



,XN 



l<z,3<N 



K 2 {x,y) = K(x,y), 

= f Si(x,y) S 1 D(x,y) 

f S 4 (2x,2y) S 4 D(2x,2y) 

\eS 4 (2x,2y) S 4 (2y,2x) 



Ki(x,y) 
K 4 (x,y) 



(A3) 

(A4a) 
(A4b) 

(A4c) 



Here det is to be interpreted as a quaternion determinant 
in the case of = 1 and 4 || . D stands for the differen- 
tiation operator, and e, Si and S4 stand for integration 
operators with kernels |^6| : 

e(x,y) = ^sgn(x- y), 

Si (a:, y) = [1 - (1 - K)eKD]~ 1 K(y, x), (A5a) 
Si{x,y) = [1- (1- K)DKe]- l K(x,y). (A5b) 

A composite operator such as eS± is defined to have a 
convoluted kernel, eS%(x, y) = J_ dz e(x, z)Si(z, y) and 
so forth, and [• ••] stands for an inverse operator. 

The probability E@ [J] of finding no eigenvalues in a set 
of intervals J is defined as 



Ep[J] 



/ dx\ ■ ■ ■ dxNVp(x\ 



,X N ). 



(A6) 
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By virtue of Eq. (A3) and the identity Yl detM n . n . oc 

det M nm , it is expressed in terms of the Fredholm deter- 
minant of the scalar or matrix kernel [ |55| : 

E 2 [J] = det(l- K 2 \j), (A7a) 
_ (A7b) 



E 1A [J] = ,/det(l- K lA \j), 



where | j represents restriction of the kernel to the inter- 
val J. 

In the following we assume that the scalar (unitary) 
kernel K is translationally invariant and symmetric (the 
latter pro pert y is respected for invariant RME's by con- 
struction (A2), but it is violated for RME's with partly 
deterministic matrix elements [p^Of), 

K(x,y) = K(y,x)=K(x-y). (A8) 

If K has these two properties, Eqs. ( |A5| ) immediately 
reduce to 

Si(x,y)=S A (x,y) = K(x-y), (A9) 

because of the relations DK = KD and eD = 1, 
and the orthogonality of the two projection operators, 
(1 — K)K = 0. We assume that liniAr^oo and the alge- 
braic manipulation that lead to the relation (A5) can be 
interchanged. In the following we let the notation K rep- 
resent its restriction to J, formerly denoted as K\ j. If J 
consists of a sing le interval [— t, t], we can simplify Ep[J] 
significantly [ml]. To do so, we introduce the resolvent 
operator R — K(l — K)~ x , and denote the associated 
kernel in the form of a matrix element, 



R(x,y) = ( x 



K 



1 - K 



(A10) 



Due to the property (A8) assumed on the kernel, it sat- 
isfies 



R{x, y) = R{y, x) = R(-x, -y). 



(All) 



Usin g this resolvent kernel, the Fredholm determinants 
(A7) are expressed as 



E 2 [-t, t] = exp | -2 J dt R(t, t) 
Ei [-t, t] = exp | - J dt[R(t, t) + R(-t, t)] 

2t 

dt[R(t, t) + R(-t,t)] 



Ei[-t,t] = - ( exp 



exp 



2/ 



dt[R(t, t) - R(-t,t)] 



(Al2a) 
(A12b) 

(A12c) 



Equation (A12a) can b e proven by taking the logarithmic 
derivative of Eq. (A7a): 



(I, nI . / 1 dK 

— lnE 2 = tr t?—t 

dt 1 ' J 1 1 - K dt 



-tr 



1 



= -2(t 



1 -K 
K 



K{\-t){-t\ + \t){t\) 
l ) ■ (A13) 



1 -K 



Eqs. (A l2b,c ) can be p roven analogously from Eqs. 
( |A4b[ ), ( |A4c|) , and ( |A7b| ). These relations are equiv- 
alent to Eqs. (6.5.19) and (10.7.5) of Ref. g because 
exp{— J Q dt[R(t, t) ± R(— t, £)]} is a Fredholm determi- 
nant of the kernel K(x, y) ± K(—x, y) (although Ref. j|] 
concerns primarily the sine kernel, the proofs of these 
equations are equally valid for any translationally invari- 
ant and symmetric kernel). 

Now we set 2t = s and denote Ep{s) = Ep[— s/2, s/2]. 
The probability Pp(s) for a pair of consecutive eigenval- 
ues to have a spacing s is clearly equal to the probabil- 
ity of finding an eigenvalue in an infinitesimal interval 
[—s/2 — e, s/2], another in [s/2, s/2 + e'] and none in be- 
tween [—s/2, s/2], divided by ee'. Thus we have 



Pp(s) = E'*(s). 



(A14) 
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FIG. 1. The decay rate n{a) = R a {po) of the LSDF. 
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FIG. 2. The LSDF P 1 (s) of the orthogonal ensemble with the kernel (|lj). 
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FIG. 3. The LSDF P 2 (s) of the unitary ensemble with the kernel 
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FIG. 4. The LSDF P 4 (s) of the symplectic ensemble with the kernel @. 
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FIG. 5. The LSDF's Pi(s) of the multifractal orthogonal ensemble (MOE) at a = 2.95 and of the Anderson Hamiltonian 
([!]) at the MIT point W = 16.4, on a lattice of size L 3 = 12 3 . Numerical data are reprinted from Fig. 1 in Ref. (23| courtesy of 
Zharekeshev. 
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FIG. 6. The LSDF's P2(s) of the multifractal unitary ensemble (MUE) at a = 3.55 and of the AH (a) under a magnetic 



field a = 1/5 at the MIT point W 
courtesy of Schweitzer. 



18.1, on a lattice of size L — 5 . Numerical data are reprinted from Fig. 1 in Ref. |24 
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FIG. 7. The LSDF's P 4 (s) of the multifractal symplectic ensemble (MSE) at a = 3.90 and of the AH (§) with spin-orbit 



coupling 6 = 7r/6 at the MIT point W = 19, on a lattice of size L 
Ref. [E5J courtesy of Kawarabayashi. 
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FIG. 8. The level number variance T, 2 (S) (divided by S) of the multifractal orthogonal ensemble at a — 3.2 and of the AH 



(^) at the MIT point W — 16.5, on a lattice of size L 6 = 10 3 . The Poisson distribution corresponds to E 2 (5')/S' : 
data are reprinted from Fig. 2 (b) in Ref. j^]] courtesy of Evangelou. 



1. Numerical 



12 



